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ABSTRACT 


Reentry of the space shuttle in the altitude region 
between 85 km and 70 km is expected to convert 1% to 10% 
of the intercepted air stream into nitric oxide. Most of 
the nitric oxide will be concentrated in a wake of small 
initial radius. The early, highly turbulent part of the 
trail will rapidly diffuse the initial high concentration 
of nitric oxide. As wake turbulence decays, atmospheric 
mixing processes become domina it and further decreases in 
contaminant concentration are slower. Atmospheric mixing 
is assumed to become dominant at a nitric oxide centerline 
concentration of 2x10 molecules/cc. Starting with these 
conditions a numerical model of eddy diffusive transport 
of emitted gases from the space shuttle wake, including 
chemical reactions between the emitted constituents and 
the ambient atmosphere, has been constructed for 75 km 
altitude. The numerical methods involve explicit solution 
of the diffusion equation and Runge-Kutta method for the 
chemical reactions. The time required to reach background 

7 

levels of nitric oxide concentration of 7x10 molecules/cc 
has been calculated. This relaxation time depends strongly 
on atmospheric conditions. Dissipation takes only a few 
hours when strong winds and a high level of turbulence 
prevail, while decay times in a quiescent atmosphere are 
measured in weeks. As the wake expands, atmospheric ozone 
interacts with nitric oxide to give molecular oxygen and 
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nitrogen dioxide. The ozone depletion in the surroundings 
of the wake results in a ♦•hole” where ozone concentration 
is less than 50% of the normal background of 10^ molecules/cc. 
For a quiescent atmosphere, this ozone hole reaches a maximum 
extent in four hours, and collapses again in less than twelve 
hours. At its maximum, the hole is one kilometer wide and 
300 meters thick, with a length of thousands of kilometers. 
Normally, strong wind shears are expected. These have been 
shown to significantly increase the rate of dispersion of 
the trail. The preliminary conclusion is that transient 
effects on atmospheric ozone are not very serious. However, 
it is recommended that the present computations be expanded 
to include typical entry trajectories, considered necessary 
because eddy diffusivity, wind shear, and atmospheric 
chemistry are all altitude dependent. The effect of en- 
hanced ionization resulting from photoionization of the 
wake should also be computed. Potential accummulation of 
nitric oxide in the mesosphere and stratosphere over time 
periods of several years at projected levels of operation 
should also be studied. 
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INTRODUCTION 


A spacecraft such as the space shuttle entering the 
Earth's atmosphere dissipates a large percentage of its 
orbital energy between 85 and 70 kilometer altitude. Some 
of this energy results in the dissociation of molecular 
oxygen and nitrogen and the formation of nitric oxide, 
which is contained in a trail of small initial radius 
e.g. 7.6 meters. The length of the trail may be as much 
as 6000 kilometers. The quantity of nitric oxide formed 
depends on several factors including spacecraft angle of 
attack, weight and surface area and can be as great as 
10% of the air swept by the vehicle. This can lead to 
initial concentrations as great as 10 14 molecules/cc at 
75 km. A discussion of this problem can be found in a 
report by Park (1972). Park estimated that shuttle flight 
operations will produce 950,000 pounds (432,000 kg) of NO 
per year. If this were to remain in a shell, 1 km thick 
at 75 km, then the background NO concentration would be 
increased by 7 x 10 molecules/cc which doubles the back- 
ground. 

The concentration of nitric oxide confined within the 
initial tube of air intercepted by the vehicle can be as 
much as 6 orders of magnitude above the average background 
density of mesospheric nitric oxide. Such large concentra- 
tions will be reduced initially by rapid mixing in the high 
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velocity near wake. However, the total time required to reach 
the background density depends on many factors incldding the 
rate of diffusion, chemical reactions and wind shear. 

In this report a numerical model of eddy-diffusive 
transport of gases in the "persistent” space shuttle wake 
is constructed. Chemical reactions between nitric oxide and 
the ambient atmosphere are included. The numerical methods 
involve explicit solution of the diffusion equation and 
the Runge-Kutta or Euler method for chemical reactions. 

Three forms of diffusive behavior are investigated, 
namely one dimensional, axisymmetr ic , and two dimensional. 

By the latter is meant a situation were vertical and hori- 
zontal diffusion coefficients are significantly different, 
but both finite. The effect of wind shear is considered 
from some of the two dimensional cases. 

The chemistry is complex and it is recognized that many 
the reactions that take place in the immediate vicinity 
of the vehicle where the air temperature is high do not 
occur in the ambient atmosphere. The next section deals 
with the chemistry and describes the reactions that have 
been included in the diffusion computations. 

CHEMICAL PROCESSES AFFECTING THE NO DISTRIBUTION 

As a nitric oxide trail expands in the mesosphere, it 
undergoes several photochemical processes. Reactions of 
importance are listed in Table I and include photodissociation, 
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reaction with atomic nitrogen and oxygen as well as the 
important interaction with ozone namely 


NO + 0 3 -» N0 2 + 0 2 


This reaction depletes atmospheric ozone by converting it 
to molecular oxygen. The effect on the ambient ozone dis- 
tribution will be large whenever nitric oxide significantly 
exceeds ozone. This is illustrated in Figure 1 where the 
ozone computed at 75 km for a nitrogen-oxygen atmosphere 
is given diurnally for different initial NO concentrations. 
In the real case the nitric oxide concentration is de- 
creasing with time at a given location due to diffusion. 
Also the reentry body has initially destroyed all ozone in 
the region where NO is produced. As will be shown, this 
’•ozone hole” expands with trail radius, as ozone reacts 
with NO to form NOg. 

ONE DIMENSIONAL DIFFUSION 

The diffusion equation is given as follows for one 
dimension : 


Sc 

FT 


- D 


5 2 c 

Sx 2 


( 1 ) 


where c is concentration, t is time, x is distance, and 
D is the molecular or eddy diffusion coefficient. In the 
above formulation it has been implicity assumed that D is 
a weak function of time and space, i.e. 
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OZONE AT 75 KM 




( 2 ) 



Now, define 

t s t'D (3) 

S 1 pi 

Since — = =• ~~r we obtain the dimensionless diffusion 
dt D ot 

equation: 


3c _ 3 C 
at ' ax 2 

The solution of the above equation as an initial value 
problem is given by Courant and Hilbert (1968) as follows: 


(4) 


+ 00 


U(x 


,t) = -~ = f 

2Vnt _J 


(?) exp(-(x-?) Z /4t) d§ 


(5) 


where i|f(?) is the distribution of concentration at time 
zero and § is the space variable. 

If the initial concentration distribution is a 
rectangle between -a and a, then equation (5) reduced to: 


a 

c(x,t) = 1 r 

c ~ 2 J 

-a 


exp(-(x-§) / 4t) d§ 


where c is the initial concentration in the interval 
o 

between -a and a . 

Equation (4) may be solved numerically, using the 
following explicit equation for marching ahead in time: 


( 6 ) 


U 


i,j+l 


= rU 


i-1, J 


+ (l-2r) U. . + rU 
> J 


i t!,j 

V 


(7 
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In the above, U is the finite difference approximation for 
concentration, the i index refers to space, and the j index 
to time. 

The quantity, r s — must be less than 0.5 in order 

(Ax ) 2 

to avoid numerical instability (Ames, 1969). 

The stability limitation can be avoided using implicit 
techniques, but such approaches are complicated and not 
easily modified. Furthermore, the stability limitation 
is not too serious. 

Let us take r = 0.4, and a diffusion coefficient, D, 

2 

of 100 m /second. 

At 7 = g- At = r (8) 

Then, the permitted time step varies with spatial grid size 
as follows: 

Ax At 7 (sec) 


1 meter 

4 x 10 

10 m 

0.4 

100 m 

40 

1 km 

1 hour 


It is readily seen from Eq. 8 that, as the diffusion 
coefficient decreases the maximum time steps become larger. 

The diffusion equation does not place any limits on 
minimum allowable time step. Thus, lower At 7 values than 
given in the table may be used if required by chemical 
kinetics considerations. 
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An example is presented in Appendix A. 

The one dimensional equation (1) can be applied to 
either horizontal or vertical diffusion. However, it must 
be pointed out that the formulation given above makes certain 
assumptions, namely: 

1. Eddy diffusion dominates 

2. Atmosphere temperature gradients are small compared 
with concentration gradients 

3. The density effect of stratification may be neglected 
for vertical distances of interest. 

Below 90 km eddy diffusion generally dominates over 
molecular diffusion. For horizontal calculations, the second 
assumption restricts us to distances small compared with these 
involved in "weather* 1 phenomena. We know that, at lower al- 
titudes, horizontal temperature gradients tend to be very 
shallow, and may reasonably assume the same to be true at 
higher altitudes, even though very little is known about 
horizontal temperature variations at altitudes above the 
cloud level. If the temperature gradients are not shallow, 
then one can expect significant winds. Generally the spatial 
and temporal history of wind patterns cannot be accurately 
determined without doing calculations of the type used in 
numerical weather forecasting. Thus, it would not be 
realistic to merely add a thermal term to the diffusion 
equation. 
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For vertical calculations, we find that the second 
and third assumptions restrict us to vertical distances of 
less than a scale height. This restriction is not serious 
for the early diffusive behavior of a trail, and is well 
justified by the resulting mathematical simplicity. The 
complete equation for vertical diffusivity is given by 
Strobel (1971) and others. 

MULTIDIMENSIONAL DIFFUSION 

One -dimensional diffusion gives a different time 
behavior than diffusion in two or three dimensions. This 
difference is of very great importance when one wishes to 
estimate residence times or lifetimes of various atmospheric 
phenomena. A comparison of the analytical formulations of 
the solution is given below (Sommerfeld, 1949) . 

The general solution is: 


u(xjt) -ymj) u( Xj t) d 
3 


(9) 


for j = 1 the integration is for one space dimension, for 
j = 2, two dimensions, and for j = 3, three dimensions. 
The kernels U for the three situations are: 


U = (4TTDt)“ 2 exp | ~^ 4 D'P"" j 
U = (4TTDt)- 1 exp { - (x - ? >/tly - ^ | 

U = (4^1) ^ exp J - 


( 10 ) 

(ID 

( 12 ) 
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The above solutions assume that the diffusion co- 


efficients are equal in all directions, an assumptions which 
is by no means always valid. 

For the sake of discussion, let us further simplify the 
above set of equations by assuming initial conditions to 
consist of simple ’’concentration poles” at the origin, and 
complete symmetry relative to the various coordinates. The 
pre-exponential factors remain unchanged, but the exponents 
become, respectively: 


-x 

4Dt 

-2x 2 

4Dt 

-3x 2 

4Dt 


(10a) 

(Ha) 

(12a) 


It is seen that the behavior of the various exponents is 
basically the same provided that the magnitude of the 
diffusion coefficient in two dimensions is one half that 
for one dimension and in three dimensions D is one third 
that for one dimension. However, the pre -exponential 
factors show a fundamentally different time dependence, 
i.e. 


1 

1 

t 

1 

v* 


for one dimension 

for two dimensions 
for three dimensions 
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One can obtain drastically different ’’lifetimes" depending 
on the number of dimensions considered. 

A comparison between one dimensional and radial dif- 
fusion is shown in Figure 2 for a diffusion coefficient 
o 

of 1000 m /sec. The initial wake radius is 10 meters. The 

12 

time required to reduce an initial concentration from 10 
8 6 

to 10 molecules/cc is 1.25x10 seconds for one dimensional 
diffusion but only 250 seconds for the radial diffusion case. 

A comparison is given between values of D differing by a 
factor of 10 for radial diffusion. The corresponding decay 
times also differ by this factor . The case where the initial 
radius is increased to 316 m is also illustrated. 

So far we have considered diffusion in a homogeneous 
medium, where the diffusivity is the same in all directions. 

The atmosphere, however, is frequently characterized by a 
situation where dif fusivities in the vertical and horizontal 
directions are significantly different. Depending on altitude 
and stratification conditions, horizontal and vertical eddy 
diffusion coefficients can differ by as much as several orders 
of magnitude. In the 75 km region, it is reasonable to as- 
sume that the ratio of horizontal to vertical eddy diffusion 
coefficients is approximately ten to one. 

It thus becomes necessary to solve the diffusion 
equation foa? inhomogeneous conditions: 

(l3) 
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It is, however, possible to examine the behavior of the 
inhomogeneous system by employing the trick of using dif- 
ferent scales in the horizontal and vertical directions. 
Re-write Eq. 13 in the form: 


where 


Sc 

at 7- 


r ,2 , .2 

= D L_£ + L_ Ls 

U 1 _ 2 + 2 - 2 

3x p 3z 


2 ^1 

6 s 4 >x 


Now, define a vertical scale, s, such that 


(13a) 


(14) 


s = pz 


(15) 


Physically, Eq. 15 means that one unit of z corresponds to 
3 units of s. Thus, if p = 100, then one unit of z, say 1 
meter, corresponds to 100 units of s, i.e. 100 cm. Mathe- 
matically, we have s = lOOz, but when dimensions are included, 
then the expression becomes s = (100 cm/m)z, i.e. the s co- 
ordinate corresponds to a finer mesh than the z coordinate. 

Equation 13 may be written in dimensionless form as: 



(16) 


where t = Dt / , and the derivatives are evaluated according 
to appropriate scales. The formal solution of Eq. 16 is 
again Eq. 11. However, the differential volume element 
considered above is a rectangle rather than the usual square. 
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In finite difference form the two-dimensional diffusion 
equation becomes 


U. . . . U. . . = r 
i,j,k+l x,j,k x 


+ r U. . , . + U. . - 2U. 


[ U i-1» J,k +U i+1 , j ,k“ 2U i, j,kj 
j+l,k " 2U i , j ,k 


where 


D, At 

^ (Ax) 2 


r = 


D 2 At 

(Ay)' 


In order to obtain the finite difference form of (16) we 


need r 


giving s’ = 5-5 1 * e * 7^ 

(Ax) 2 (Ay) 2 “ 


Again, the direction with the small diffusion co- 
efficient requires a finer mesh. 

It appears from the above disucssion that one may 
examine concentration profiles in both directions, while 
only using a fine mesh in the direction where diffusion is 
slow, and gradients steep. If, however, one is not very 
interested in the detailed development of the vertical 
profile, then it may be just as well to use equation (17) 
with Ax = Ay, in which case 


r = (^) r 
y 'd , 1 x 


While proper treatment of diffusion in the atmosphere 
frequently requires use of a two-dimensional calculation, 
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it is interesting to examine which conditions permit a 
simplified treatment. 

To facilitate visualization of when to make the 
transition from one -dimensional treatment to two-dimensional 
treatment, re-write equation (11) in the form: 


where 


U = AB 


_1 

A = (4nD^t)~ 2 exp 


_i 

B = (4nD 2 t)"" 2 exp 


-(x-;) 2 

4D ± t 

-Cy-T]) 2 

4D 2 t 


(19a) 


(19b) 

(19c) 


Now lim (B) = 1. Thus, in dimensionless form the criterion 
D2t-»o 

for using the one -dimensional formulation is: 


D 2 t 


(Ay) 


« l 


( 20 ) 


where Ay is the initial vertical width of the uniform con- 
centration regime. 


DIFFUSION WITH CHEMICAL REACTION 

One dimensional diffusion with chemical reaction is 
described by the equation: 






where is the rate of loss of specie "i” . 


N 

R . = £ k . c . c. 
i J= i j i k 


( 21 ) 


( 22 ) 
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is the net loss of specie "i” by the relevant chemical 
reactions j with species "k" ; j and k go from a minimum of 
one to a maximum of N. 

The atmosphere is treated as being isothermal, so that 
the temperature only comes in at the initial evaluation of 
the rate coefficient. 

One may break up Eq. (21) into two components, i.e. 


where 


and 


3c l . i^iP * SC iC 

st 7 at • at 7 


dC 


iD 


st‘ 


D i 


a 2 c. 


Sx 


9c iC N 

= -£ k, c, c, 
at j=i 3 1 b 


(23) 


The space field can be swept for chemical reaction 
during a time step, At, followed by a sweep for diffusive 
relaxation for that same time step. In this way it is pos- 
sible to use the explicit marching formula for diffusive 
relaxation, and Runge-Kutta for the solution of a set of 
coupled nonlinear first order equations. 

To check the validity of the above procedure a compari- 
son was run between the method described above, and a direct 
solution of Eq. (21) as written. No difference could be 
seen in the results. The Euler method was used for the 
chemical component in both cases. 
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Computations were carried out for one dimensional 
diffusion with kinetics, radial diffusion with kinetics, 
and two dimensional diffusion with kinetics. The rationale 
behind the computational methods has already been discussed 
and will not be repeated here. Solutions to the two 
dimensional diffusion equation with chemistry are pre- 
sented for illustrative purposes both as three dimensional 
plots, and as contour plots. In both situations a series 
of "frames" developed in time are shown. Figure 3 shows 
a three dimensional projection of the three chemical 
species. The horizontal diffusion coefficient is 
l(j cm /sec, and the vertical coefficient is 10 cm /sec. 
Figure 4 shows contour plots for a similar situation. 
APPLICATION TO THE SPACE SHUTTLE 

A space shuttle reentry schematic is shown in Figure 5. 
Three basic regimes may be distinguished. First is the 
"near" wake mixing regime where eddy diffusivity is deter- 
mined by turbulence generated by shear due to the high reentry 
velocity. This regime has been treated by Park (1972) who 
used a model developed by Lin and Hays (1965) . The Lin-Hays 
model is quasi-one dimensional, and primarily designed as 
a method to treat wake chemistry without regard to radial 
wake gradients which are assumed to be zero. 

In using this model Park assumes that the initial wake 
radius for the shuttle is 7.6 m and that its boundary varies 
according to the formula 


16 






18 






19 



( 24 ) 



r 

o 


1/3 


0.7 (S_) 
v r ' 
o 


where r fo is the wake radius at any point, r Q the initial 
wake radius and 


s' = s Q / + x (25) 

is the distance from a hypothetical apex of the wake, x is 

the distance from the section where r^ = r Q . s q / is the 

distance from the hypothetical apex to the section where 

r, = r . For r = 7.6 meters, s ' = 21.6 meters. Also, 
boo ’ o ’ 

define a quantity i/ as the distance from the leading edge 

of the vehicle to the section where r. = r . Reference to 

b o 

Figure 9 in Park (1972) , reproduced in the Appendix as 
Figure A-4, shows that i/ is approximately 91.5 meters. 

Park also defines a quantity, s as the distance from the 
leading edge of the vehicle. We have 

s = i/ + x (26) 

Consider a wake radius, r of 23 meters, giving 
r b ° 

— = 3. Calculating back from Eq. 24 one obtains s' = 600 
1 o 

meters. From Equations 25 and 26 we then deduce that s, 
the distance from the leading edge, is 670 meters (2200 ft.). 

Reference to Park’s Figure 10, which shows the amount 
of air converted to nitric oxide as a function of distance 
from the leading edge of the vehicle at 70 km indicates 
that 670 meters is the section at which the ratio of nitric 
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oxide mass produces to mass of air swept is at a maximum. 

The ratio varies between approximately 2.5% and 9% depending 
on vehicle entry velocity. The mole fraction variation is 
approximately the same. According to Goldburg (1965) 

610 m is approximately the point where transition from a 
laminar to a turbulent wake occurs. Thus, 670 m is some 
60 m into the turbulent wake regime. The value at 670 m 
is 38, where D = 2 r^. 

Reference to Bailey (1966) , Li (1963) and Slattery 
and Clay (1962) shows that at the section considered above, 
the wake centerline velocity is still some 10% of vehicle 
velocity, i.e. in the order of 610 m/sec. 

However, the wake centerline velocity is down to 1% 

of vehicle velocity, i.e. 61 m/sec at - 420. Such a 

value is comparable to wind speeds at 75 km altitude 

(Smith et al., 1964). The above value corresponds to 

4 

s = 7,400 m giving r fe = 52 meters. At s = 3x10 m the wake 
velocity would decrease to 9.1 m/sec (Park, 1972) which is 
the same order as the velocity of natural turbulence. At 
this point the wake radius, r^, would be 132 meters. 

In the ’'near** wake mixing situation discussed above 
eddy diffusivity is determined by turbulence generated by 
shear due to the high Teehtry velocity. This regime may 
be treated using a time dependent eddy diffusivity whose 
form can be deduced by matching to the empirical wake growth 
law or in the manner of Park (1972) . The use of a time 
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dependent eddy diffusivity combined with chemistry is 
actually superior, as it takes account of radial gradients 
existing in the wake. These gradients are assumed to be 
zero in the Lin and Hayes (1965) model used by Park. The 
time dependent eddy diffusivity approach is an interesting 
possibility for future study. 

For the present we have assumed that the Park-Teare 
calculation is valid up to a certain point in the wake, 
after which the dominant process becomes atmospheric mixing 
with constant eddy diffusivity. The proper point of transi- 
tion from wake generated turbulence to atmospheric turbulence 
is not constant, but depends very much upon atmospheric 
conditions. The more turbulent the atmosphere, the sooner 
atmospheric mixing becomes dominant, and vice versa. 

Referring once again to Figure 5 the second basic 
regime may be termed the atmospheric mixing regime. Chemical 
reactions are important, but affect eddy diffusive processes 
primarilty through modification of the boundary conditions 
for the diffusive process. 

In the third regime, concentration gradients have 
become very shallow, and the rate of diffusion slow. Chemical 
rate processes now play the dominant role. This regime pro- 
vides proper input parameters for more sophisticated chemical 
reaction models of the atmosphere. 

Figure 6 shows the centerline nitric oxide concentra- 
tion for a reentering shuttle. Initial conditions for the 
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NEAR WAKE REGIME 


ELLIPTICAL DIFFUSION 
D hor = 10 5 cm 2 /sec 

4 2 

D yER = 10 cm /sec 

k = 0.194 * 10’ 14 /sec 


ALTITUDE = 75km 

R -90 METERS 



CHEMISTRY DOMINATES 





atmospheric calculation were taken as follows: nitic oxide 

concentration is taken as 200 x 10^ molecules/cc, and 

initial radius as 90 meters. The ratio of horizontal to 

vertical eddy diffusivity is taken to be ten, the diffusivity 

5 2 *~1 2 

coefficients being 10 cm sec (10 m /sec) horizontal and 
lO^cm^sec -1 (1.0 m^/sec) vertical. These are conservative 
estimates, and cases with other values of eddy diffusivity 
will be discussed later. Extrapolations for one dimensional 
and radial diffusion are shown, and transition to the chemi- 
cally dominant regime is seen to occur after approximately 
one day. Figure 7 shows development of concentration pro- 
files with time in both the vertical and horizontal directions. 
In this case we started with a situation where initial ozone 
concentration grows monotonically with time. A more inter- 
esting quantity in this case is the "ozone hole”. Figure 8 
shows development of the half width of the point of 50% 
ozone depletion, i.e. inside the radius shown ozone concen- 
tration is less than 50% of the natural ambient value. For 
the case considered in Figures 6 and 7 we see in Figure 8 
that the maximum horizontal spread of the ozone hole occurs 
in approximately 3 hours, and amounts to a 400 meter half 
width for the 50% point. After approximately 5 hours, the 
gap begins to close, and after 10 hours the 50% point has 
moved to the center of the trail. At 20 hours, centerline 
ozone concentration is 85% of the ambient value, whereas 
centerline nitric oxide concentration is still 20 times as 
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high as the natural background. The vertical extent of the 
ozone depletion is less, since the vertical diffusion co- 
efficient has been taken to be only 10% of the horizontal 
value. Otherwise, vertical behavior is quite similar to 
horizontal behavior. 

Thus far, we have only discussed what may be considered 
a ’’night time" situation, i.e. nitric oxide reacts with ozone 
to give nitrogen dioxide whose further chemical activity is 
not considered. However, much of the alarm about nitric 
oxide pollution has referred to the "daytime" situation in 
which nitrogen dioxide is photodissociated to nitric oxide, 
thus creating a catalytic cycle in which ozone is destroyed 
without any corresponding reduction in odd nitrogen. The 
purely photochemical case is serious indeed. However, 

Figure 8 shows that in the diffusive situation the recycling 
of nitric oxide has a relatively minor effect. The peak 
half width of the ozone hole is increased by only 5%, while 
the collapse time goes from 32,000 seconds to 48,000 seconds, 
i.e. from 9 hours to 13 hours. Figure 9 shows nitric oxide 
decay for both night and day situations. The daytime nitric 
oxide does decay more slowly as expected, but the difference 
between daytime and nighttime cases only becomes significant 
for long times. It might be noted that, for the daytime 
situation, nitric oxide concentration is govered only by 
diffusion, since any nitric oxide which reacts with ozone 
is almost immediately recycled because of the rapid 
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photodissociation rate. In the situation considered, it 
was assumed that ozone concentration is zero in the space 
originally occupied by the nitric oxide. Thus, ozone con- 
centration increases monotonically to some asymptotic value. 

At the other limit, where background ozone concentration 
extends into the core of the initial high nitric oxide 
concentration regime, a rapid chemical reaction occurs, 
bringing ozone in the core region to a very low concentration, 
i.e. a situation very similar to the one considered. We see 
from Figure 10 that the initial increase in centerline ozone 
concentration is very rapid, a behavior attributable to 
initially steep concentration gradients and small effect 
of chemistry. In half an hour a plateau region is reached 
which lasts for approximately four hours. The plateau is 
the regime of most active chemical reaction, as is also 
evident from the nitrogen dioxide plot which reaches a 
maximum in this time period. A look at the rate constant 
for the reaction: 

NO +■ 0 3 -> N0 2 + 0 2 

explains the behavior obtained. The rate constant for the 

-4 

above reaction is 0.2x10 seconds when concentrations are 
measured in terms of lO 1 ^ molecules per cubic centimeter. 

When the plateau is reached, the ozone concentration is 
0.3 in the above units, while nitric oxide concentration 
is 10. Thus, for the above chemical reaction: 
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Ozone recovery (centerline) 




C^° 2 1 = (0.2x10 4 sec -1 ) (10) (0.3) = 0. 6xl0“ 4 sec _1 

4 

giving a characteristic time of approximately 2x10 seconds, 
which is indeed seen to be characteristic for the nitrogen 
dioxide buildup. As the ozone buildup leaves the plateau 

4 

region at 16x10 seconds, the nitric oxide concentration 

4 

has already decreased to 4.5, and at 34x10 seconds the 

nitric oxide concentration is down by a factor of ten to 

1.0, while ozone has only increased by a factor of two, 

thus decreasing the nitrogen dioxide production rate to 
-4 -1 

0.1x10 sec . However, the rate of outward diffusion has 
only increased, since centerline nitrogen dioxide concentra- 
tion is greater than before. Thus, the nitrogen dioxide 
concentration begins to drop, and continues to do so, since 
the rate of production continues to decrease. 

Reference to Figure 10 shows that nitrogen dioxide 
concentration is always below 0.7xl0 10 molecules/cc. Thus, 
one would expect the photodissociation back to nitric oxide 
to become important only after nitric oxide concentration 
becomes comparable. References to Figure 8, 9 and 10 
however, indicate noticeable difference between night and 

4 

day situations as early as 6x10 seconds, which is the point 
at which maximum nitrogen dioxide concentration is approached. 
Furthermore, it is found that the daytime nitric oxide con- 
centration is the sum of nighttime nitric oxide and nitrogen 
dioxide, at least for time in excess of half an hour. A 
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steady state analysis of the daytime chemistry shows that 
nitrogen dioxide concentration is approximately 0.2% of the 
nitric oxide value. This is quite different from the night- 
time situation where nitrogen dioxide concentration actually 
exceeds nitric oxide concentration for times longer than 
15 hours. 

The cases discussed above all have diffusion coefficients 
of 10 5 cm 2 /sec horizontal, and 10 4 cm 2 /sec vertical. These 
estimates are conservative, representing reasonable lower 
limits to diffusion coefficients which may be expected at 
75 km altitude. However, it is essential to recognize that 
knowledge of atmospheric diffusion coefficients is very 
sparse, and that large variations in eddy diffusivity are 
quite probable. 

For the daytime situation, where the effect of chemistry 
on nitric oxide concentration is slight, one may expect that 
the time required to reach background concentration levels is 
inversely proportional to the diffusion coefficient. One 
obtains approximately the same behavior for nighttime 
centerline nitric oxide concentration, although the propor- 
tionality is no longer quite accurate. For centerline ozone 
concentration such simple proportionality generally does not 
hold, and a more complete analysis is needed. 

To study the effect of eddy diffusivity, we ran three 
cases in which the initial conditions were identical, the 
number of time steps were the same, but the eddy diffusion 
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coefficients were varied. The situation of homogeneous 
radial diffusion was chosen for this comparison. The 
initial radius was taken to be 20 meters, and the initial 
nitric oxide concentration was 10^ 4 molecules/cc. Results 


of the above 

comparison 

are given in 

the following 

table : 

o 

D,m /sec 

t 7 (sec) 

Dt'/r o 2 

[NO] 

[0 3 ] 

100 

2xl0 3 

500 

6xl0 10 

0 . 8xl0 10 

10 

2xl0 4 

500 

5xl0 10 

0.21xl0 10 

1.0 

2xl0 5 

500 

5xl0 10 

0.2xl0 8 


where r is initial radius and t' is time after 20,000 
o 

intervals. It is seen that for constant values of the 

n 

quantity (Dt 7 /r Q ) the centerline nitric oxide concentration 
is approximately the same, but that centerline ozone con- 
centration decreases for decreasing diffusion coefficient. 
The reason for the above behavior is that the slower the 
diffusion, the longer the chemical destruction of ozone can 
proceed before nitric oxide concentration drops to a low 
value. A comparison of concentration -time histories is 
shown in Figure 11. The eddy diffusion coefficient is 
varied. 

Figure 12 shows a comparison of one dimensional two 
dimensional, and radial diffusive behavior. The one dimen- 
sional centerline concentration decay approaches at 2 
asymptote, and the radial decay has a t -4 asymptote. For 
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TIME, SECONDS 

Dependence of centerline concentration on diffusivity 








the two dimensional case, we took the ratio of vertical 
to horizontal eddy diffusivity to be 0.1. The graphs on 
Figure 12 show that the two dimensional case is initially 
close to one dimensional behavior, then approaches a radial 
type time decay for long times. Reference to equations 19 
and 20 shows that the above behavior is precisely what one 
would expect from theory. There, it was found that one 
dimensional behavior would be observed when: 

4D 9 t 7 

p = _2 << x 

V 

where Dg is vertical eddy diffusion coefficient. From 
Figure 12 we see that at t' = 1,000 seconds deviation from one 
dimensional behavior is 17% and at that point the quantity p 
is equal to 0.4. A more general examination of the dif- 
ferences between the two curves shows that, very roughly, 
the deviation from one dimensional behavior is equal to 
(p/n). Radial time decay is reached when P > 1. However, 
the effective diffusion coefficient is some mean of vertical 
and horizontal coefficients. For the particular case chosen 
above, the best mean is simply the arithmetic average of 
the coefficients in the two directions, but such simplicity 
is seen to be fortuitous when other ratios of vertical to 
horizontal diffusivity are considered. 

THE EFFECTS OF WINDS 

Thus far, winds have been considered implicitly, as 
generators of the background turbulence which gives rise 
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to the eddy diffusivity. Winds also affect the persistent 
wake in two other ways, namely advection and distortion. 

The former is, of course, due to a steadily blowing wind, 
while the latter is caused by wind shear. 

First, consider the effect of advection. Berman and 
Goldburg (1971) quoted an annual mean zonal wind velocity 
of 10 meters/second for the stratosphere. The work of 
Smith et. al., (1966) suggests higher wind speeds at higher 
altitudes. Taking 10 m/sec as a reasonable lower limit, 
one finds that in 10,000 seconds (approximately 3 hours) 
a North-South trail will have been wafted some 100 km in 
the direction of the prevailing wind. A trail originally 
deposited in the East -West direction will have been moved 
along its own axis, giving little indication of wind effect. 
The above argument assumes a purely zonal wind, while a 
wind having both zonal and meridional components is prob- 
able more realistic. In that case, the corresponding 
effects are determined by the proper components of the 
wind vector relative to the wake axis. 

Extrapolation of higher altitude wind shear data of 
Bedinger and Constant inides (1969) to 75 km indicates an 
average vertical shear of 7 meters/second/kilometer, or 
0.7 meters/second/100 meters. Thus, two horizontal layers, 
separated in altitude by 100 meters, would move at an? average 
rate of 0.7 meters/second relative to each other. In three 
hours, the displacement would therefore be as large as 7 km. 
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If the wind shear is perpendicular to the trail axis, then 
the trail cross section will be stretched out into a thin 
oval shape, decreasing the thickness of the ozone hole, and 
increasing the rate of dissipation of nitric oxide. On the 
other hand, a trail whose axis is parallel to the axis of 
shear will only suffer a slight distortion in the vertical 
concentration profile. 

Figure 13 shows the effect of wind shear perpendicular 
to the wake axis. Compare Figure 13 with the vertical section 
of Figure 7. Note that the effect of shear is to slant the 
concentration profiles as gas from the periphery is moved 
closer to the center. Figure 14 gives a comparison of 
centerline nitric oxide concentration with, and without 
wind.-sbe?r. Nighttime chemical conditions were chosen. It 
is readily seen that wind shear plays quite an important 
role, and should be studied in greater detail. 

EFFECT OF THE BOW SHOCK 

Park (1972) considered the nitric oxide produced in the 
stream tube which becomes the wake of the vehicle, and the 
present report has been concerned with the transient effect of 
the wake nitric oxide. However, there-is additional nitric 
oxide produced by the extended bow shock. To properly 
estimate the amount thus produced one needs to know the 
exact shape of the bow shock, as well as the flow 
characteristics in the stream tubes behind it. The above 
requires a detailed aerodynamic calculation not available 
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Effect of vertical shear 




to the author. For a rough estimate it was assumed that the 

13 

bow shock produces a nitric oxide concentration of 10 

molecules/cc (i.e. 1%) over an initial radius of 400 meters. 

2 

The horizontal diffusivity was taken to be 100 m /sec, 

2 

and the vertical diffusivity was 10 m /sec. It was found, 
that in the absence of shear the ozone hole after 12 hours 
is some 3.5 km wide, and centerline ozone concentration is 
0.3% of background. However, introducing a perpendicular 
shear of 1.25 m/sec/km made the ozone hole close after 
4 1/2 shours. 

The above calculation is only a rough estimate, but it 
shows the importance of a detailed calculation of the bow 
shock production of nitric oxide. 

ENHANCEMENT OF THE LOWER IONOSPHERE BY NITRIC OXIDE INJECTION 

The nitric oxide and nitrogen dioxide resulting from 

the reentering body, can be ionized by direct solar Lyman 

alpha radiation during the day and by scattered Lyman alpha 

12 —3 

radiation at night. A starting NO density of 10 cm at 

6 3 

75 km will give rise to 2x10 electrons/cm for a zenith 
angle of 7.5°. This is more than three orders of magnitude 
above the normal background. The majority of the primary 
N0 + and N0 2 + ions formed dissociatively recombine with free 
electrons. However, ionic reactions do occur which modify 
the primary species. For instance NC> 2 is transformed into 
N0 + by the reactions 

N0 2 + + NO -* NO + + N0 2 
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The ion N0 + can react with NO, C0 2 or N 2 to form a cluster 
ion which in turn is transferred into a hydrated ion of the 
form H q 0 + *(H o 0) . Such hydrated ions may further cluster to 
form an aerosol under proper conditions of humidity, ; pressure 
and temperature. The dissociative recombination rate of 
hydrated ions is as much as a factor of 10 larger than that 

-j- "|* 

of NO or N0 2 . In general the electron density is great 
enough that NO and N0 2 will directly recombine but as 
the electron density decreases the percentage of hydrated 
ions will increase. 

REACTIONS BETWEEN HEAT SHIELD ABLATED MATERIAL AND THE 
ATMOSPHERE 

It is understood that if an ablative heat shield is 
used one possible product of ablatement is Si. This 
material can form long-lived ions. These ions normally 
recombine with electrons by the process of radiative 
recombination 

Si + + e Si + hv 

-12 3 —1 

and have a rate coefficient of approximately 1x10 cm sec 
which yields a long lifetime. In addition to photoionization 
Si can be created by charge transfer with other molecular 
ions by reactions such as 

N0 + + Si -* Si + + N0 2 

and 

N0 2 + + Si Si + + N0 2 
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Molecular ions involving Si can be formed by associa- 
tion 

Si + + XO -> SiO + + X 

where XO is a neutral such as Og. The ion SiO + would then 

-7 3 -1 

recombine at a rate of the order of 4x10 cm sec , which 
is orders of magnitude faster than radiative recombination. 
However, the reaction 

SiO + + O -> Si + + 0 2 + 0.6 ev 

also occurs with a rate coefficient of 2xl0 - *"® cm** sec - * - 
at 300°K. 

CONCLUSIONS AND RECOMMENDATIONS FOR FURTHER STUDY 

Reentry of the space shuttle raises the concentration 
of nitric oxide by five orders of magnitude in a wake whose 
initial radius is approximately 8 meters. A relatively per- 
sistent wake is formed. A numerical model has been constructed 
of the eddy diffusive dispersion of the wake. Chemical re- 
actions between wake nitric oxide and the ambient atmosphere 
have been included, particularly the reaction between nitric 
oxide and ozone which results in the formation of nitrogen 
dioxide . 

The rate of dispersion of the wake depends on the 

magnitude of the eddy diffusion coefficient and wind 

structure. For radial diffusion and an eddy diffusion 

2 

coefficient of 1100 m /sec the centerline nitric oxide 
concentration is reduced 4 orders of magnitude in 250 
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seconds. While such high diffusivity may be expected in the 

extremely turbulent early wake regime, the natural atmospheric 

diffusivity , is considerally ! less. Thu§, Johnson .and a W i Ik ins 

(1966) estimate the upper limit of vertical eddy diffusivity 

5 2 2 

to be approximately 8x10 cm /sec, i.e. 80 m /sec. Since 

atmospheric eddy diffusivity may be expected to vary by 

several orders of magnitude we have chosen to consider 

2 

a vertical diffusivity of 1 m /sec as a reasonable lower 

limit. It must be noted, however, that molecular diffusivity 

2 

at 70 km altitude is approximately 0.5 m /sec. Thus, we 

may be reasonably confident that our calculation is, indeed, 

2 

conservative. In the case of 1 m /sec the time to reach 
background can be weeks. .Normally strong wind shears 
are expected. When present these can significantly in- 
crease the rate of dispersion of the trail. 

As the wake expands atmospheric ozone is converted to 
molecular oxygen and nitrogen dioxide by nitric oxide. The 
ozone depletion results in an ’’ozone hole” in the surroundings 

of the wake. For a vertical eddy diffusion coefficient of 
2 

1 m /sec and no wind this ozone hole reaches its maximum size 
in approximately four hours and collapses again in less than 
12 hours. At its maximum the hole is approximately one 
kilometer wide and 300 meters thick with a length of 
thousands of kilometers. 

The preliminary conclusion is that the transient effects 
of a reentering shuttle on atmospheric ozone are not very 
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severe, and do not last much longer than a day. However, 
large numbers of shuttle flights do pose a potential ac- 
cummulation of nitric oxide which may seriously affect ozone. 

It must be recognized that all computations depend on 
a number of parameters which are ill -defined. A thorough 
knowledge of the value of the eddy diffusion coefficient 
as a function of position and time is essential. The same 
holds true for wind profiles. The lifetimes of the trail 
is a sensitive function of these parameters and there are 
few experimental values. Other quantities such as the 
ambient distribution of nitric oxide and the ionospheric 
properties should also be studied in order to better under- 
stand the physical and chemical properties of the atmospheric 
entry altitude regime. 

It is recommended that the present computations be 
expanded to include typical entry trajectories, considered 
necessary because eddy diffusivity, wind shear and atmo- 
spheric chemistry are all altitude dependent. Thus, 
significant differences in contaminant behavior may be ex- 
pected. The effect of enhanced ionization resulting from 
photoionization of the wake should be computed. Nitric 
oxide production in the extended bow shock should be con- 
sidered in detail. 

A series of nitric oxide release experiments should 
be conducted in the mesosphere under different atmospheric 
conditions . 
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APPENDIX A 


A one dimensional problem was solved as followed: A 

length L was considered in which an initial rectangular 
concentration profile of width l/ = 0.1L was located in 
the center. A history of the centerline concentration as 
a function of the number of time steps is given in Figure Al. 

Figure A2 shows growth of the width at half the center - 
line concentration where X is the half width. The asymmetry 
for long times is believed to be largely due to the way the 
boundary conditions were chosen, i.e. the concentration at 
the boundary of the intervals was taken as a fraction of 
the concentration at a fixed point close to the boundary. 

Now, consider application of the calculation to the 
earth's atmosphere. 

From Eq. 8 in the main text we may conclude that 


/ 

t 


nAt ' 


PAt _ (Ax ) 2 

D “ D 


(A-l) 


where n is the number of computation steps. 

At 75 km altitude, Zimmerman and Champion (1963) give 

an effective diffusion coefficient in the range of approxi- 

2 

mately 300-2000 m /sec. Johnson and Wilkins (1966) deduce 
2 

100 m /sec. To be conservative, let us choose the lower 

2 

value of D, i.e. 100 m sec. 

The numerical solution discussed above shows that it 
takes 100,000 time steps for the centerline concentration 
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250 
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n , THOUSANDS 

Half width versus number of time steps 





to drop from 1000 units to 18 units. Reference to Eq. 4 
in the main text shows that concentration appears to the 
same power on both sides of the equation. Therefore, c may 
be multiplied by any constant without altering either the 
differential equation or its solution. Thus the initial 
1000 units may be set to any initial concentration c Q , or 
mole fraction, M Q . 

In the numerical computation, the total lateral extent, 
L, was broken up into 400 intervals, i.e. Ax = the 

ratio, r, was taken as 0.4. 

2 L 2 

Then, using Eq. A-l: (Ax) = giving 

16x10 

2 

t' = ( n \ rh n 

X ''10,000' MOD' 

Equation A-2 may be used to estimate the actual rate of 
diffusion in the atmosphere. Consider the time, t*, re- 
quired for the centerline concentration to diminish from 
1000 units down to 18 units. We know from the computation 
that this required 100,000 time steps. Thus we have 



It is now a simple matter to estimate t" for various 
initial widths of the high concentration peak. D is taken 
as 100 m^/sec, and %" is given for several i/ and L values, 
i/ is the width of the original concentration peak. 


(A-2) 


(A-3) 
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L t "(jjec) t" (D=100) t" (D=300) t"(D=1000) 

4m 40 m 4 4 sec 1.3 sec 

40 m 400 m 400 6.7 min 2.3 min 

1/2 mile 8x10 m 16x10^ 44.5 hours 15 hours 

* D = m /sec 

A way to check the validity of the above computation is 
to compare the predictions of Figure A**2 with cloud spreading 
data measured by Zimmerman and Champion. These authors show 
that a contaminant cloud at 75 km requires approximately 
2000 seconds to grow from 800 meters diameter to 8000 meters 
diameter. Figure A-3 shows that the width at half maximum 
grows by a factor of 10 in 13,500 time steps. Equation A-2 
may now be used to calculate an effective diffusion co- 
efficient. L - 10 i/, n - 13,500, t' = 2 x 10 3 sec. Thus 


.4 sec 
40 sec 
4.5 hours 


2xl0 3 


a. 35) (■■ 8 ^° 3]2 ) 


1 

D 


2x10 x 40 


(8x10°) (8x10; (1.35) 


It follows that D ~ 1080 m 2 /sec. 


Now the effective diffusion values quoted by Zimmerman 

and Champion for 75 km range between 300 and 2000 meters 2 /sec, 

2 

so that the 1100 m /sec calculated above is probably more 

o 

reasonable than the 100 m /sec used previously. Nontheless 
it must be remembered that atmospheric eddy diffusion can 
easily vary by one, or even two, orders of magnitude depending 
on the activity of the atmosphere. 
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WAKE BOUNDARY 
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used in computer calculation, from Park (1972) 


TABLE I 

NEUTRAL REACTIONS 


REACTION 


0 2 + hv 

-> 

0 + 0 

°3 + hv 


o + o.jC 1 

0 + 0 2 + M 

-> 

0 3 + M 

b + o 3 


°2 + °2< 

NO + hv 


N + 0 

N + 0 

— > 

NO 

N 2 + hv 

-» 

N + N 

NO + 0 3 

-* 

N0 2 + 0 2 

N + N0 2 

-> 

2 NO 

0 + N0 2 

-> 

NO + 0 2 

N0 2 + hv 


NO + 0 

NO + NO + 0 2 

-» 

N0 2 + NO, 

NO + 0 + M 


N0 2 + M 

NO + 0 

-» 

N0 2 

NO + 0 2 

-> 

N0 2 + 0 

N0 2 + no 2 

-» 

NO + NO 


S' 




RATE 


8( -35) EXP(445/T) 

5 . 6( -11) EXP( -2850 /T) 

J NO 

1(-17) 

J N 

1 ( -12) EXP( -1200/T) 

6 ( - 12 ) 

1 . 66(-ll)EXP( -300 /T) 
J NO z 

6 . 61 ( -39) EXP( 526/T) 

3 . 63 ( -33) EXP(930/T) 

4. 07 (-12) 

1 . 58 ( -12) EXP( -22900/T) 
6 . 61 ( -12) EXP( -13500/T) 
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